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Abstract. We study the long-term evolution of massive black hole binaries (MBHBs) at the 
centers of galaxies using detailed full three-body scattering experiments. Stars, drawn from 
a distribution unbound to the binary, are ejected by the gravitational slingshot. We quantify 
the effect of secondary slingshots - stars returning on small impact parameter orbits to 
have a second super-elastic scattering with the MBHB - on binary separation. Even in the 
absence of two-body relaxation or gas dynamical processes, very unequal mass binaries 
of mass M = 10 7 M Q can shrink to the gravitational wave emission regime in less than a 
Hubble time, and are therefore a target for the planned Laser Interferometer Space Antenna 
(LISA). Three-body interactions create a subpopulation of hypervelocity stars on nearly 
radial, corotating orbits, with a spatial distribution that is initially highly flattened in the 
inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. 
The mass ejected is > 0.7 times the binary reduced mass, and most of the stars are ejected 
in an initial burst lasting much less than a bulge crossing time. 
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1. Introduction 

It is now widely accepted that the formation 
and evolution of galaxies and massive black 
holes (MBHs) are strongly linked: MBHs are 
ubiquitous in the nuclei of nearby galaxies, and 
a tight correlation is observed between hole 
mass and the stellar mass of the surround- 
ing spheroid or bulge (e.g. Magorrian et al. 
1998; Gebhardt et al. 2000; Ferrarese & Merritt 
2000; Haring & Rix 2004). If MBHs were also 
common in the past (as implied by the no- 
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tion that many distant galaxies harbor active 
nuclei for a short period of their life), and if 
their host galaxies experience multiple merg- 
ers during their lifetime, as dictated by popular 
cold dark matter (CDM) hierarchical cosmolo- 
gies, then close MBH binaries (MBHBs) will 
inevitably form in large numbers during cos- 
mic history. A MBHB model for the observed 
bending and apparent precession of radio jets 
from active galactic nuclei was first proposed 
by Begelman, Blandford, & Rees (1980). The 
coalescence of two spinning black holes in a 
radio galaxy may cause a sudden reorienta- 
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tion of the jet direction, perhaps leading to the 
so-called "winged" or "X-type" radio sources 
(Merritt & Ekers 2002). Recently, observations 
with the Chandra satellite have revealed two 
active MBHs in the nucleus of NGC 6240 
(Komossa et al. 2003), and a MBHB is inferred 
in the radio core of 3C 66B (Sudou et al. 2003). 

BH pairs that are able to coalesce in less 
than a Hubble time will give origin to the 
loudest gravitational wave (GW) events in 
the universe. In particular, a low-frequency 
space interferometer like the planned Laser 
Interferometer Space Antenna (LISA) is ex- 
pected to have the sensitivity to detect nearly 
all MBHBs in the mass range 10 4 - 10 7 M 
that happen to merge at any redshift during the 
mission operation phase (Sesana et al. 2005). 
The coalescence rate of such "LISA MBHBs" 
depends, however, on the efficiency with which 
stellar and gas dynamical processes can drive 
wide pairs to the GW emission stage. 

Following the merger of two halo+MBH 
systems of comparable mass ("major merg- 
ers"), it is understood that dynamical friction 
will drag in the satellite halo (and its MBH) to- 
ward the center of the more massive progenitor 
(see, e.g., the recent numerical simulations by 
Kazantzidis et al. 2005): this will lead to the 
formation of a bound MBH binary in the vio- 
lently relaxed core of the newly merged stel- 
lar system. As the binary separation decays, 
the effectiveness of dynamical friction slowly 
declines because distant stars perturb the bi- 
nary's center of mass but not its semi-major 
axis (Begelman, Blandford, & Rees 1980). The 
bound pair then hardens by capturing stars 
passing in its immediate vicinity and ejecting 
them at much higher velocities (gravitational 
slingshot). It is this phase that is considered the 
bottleneck of a MBHB's path to coalescence, 
as there is a finite supply of stars on intersect- 
ing orbits and the binary may "hung up" before 
the back-reaction from GW emission becomes 
important. This has become known as the "fi- 
nal parsec problem" (Milosavljevic & Merrit 
2003). 

Recently, in collaboration with Piero 
Madau and Monica Colpi, the group in Como 
started a study of MBHB dynamics from two 
different perspectives, namely, the role of gas 
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in the dynamical friction regime (the 100-to-l 
pc decay), and the role of 3-body encounters 
with bulge stars. The interested reader can find 
exahustive descriptions of computational tools, 
and detailed discussions of the results in Dotti, 
Colpi & Haardt (2006), and in Sesana, Haardt 
& Madau (2006). We have tried to answer a 
number of questions, i.e., 

i) Is stellar slingshot able to drive MBHBs 
to the GW regime in an Hubble time? What is 
the role of different mass ratios and eccentrici- 
ties in setting the orbital decay time scale? 

ii) What are the cinematical properties of 
the stellar population after the interaction with 
the MBHB? Hyper velocity stars can be an ob- 
servable large scale signature of a past MBHB 
merging. 

iii) How do eccentric orbits evolve in a 
gaseous circum nuclear disk? Do they become 
circular? This issue may be relevant in estab- 
lishing the initial conditions for the braking 
of the binary due to the slingshot mechanism 
and/or gaseous gravitational torque. 

iv) During the sinking process, do the BHs 
collect substantial amounts of gas? This is a 
query related to the potential activity of a MBH 
during a merger and its detectability across the 
entire dynamical evolution. This issue is re- 
lated to the search of EM counterparts to GW 
signals. 

In the following, I briefly review our main 
results concerning the hardening in a time 
evolving stellar background, i.e., I will discuss 
our answers to questions i) and ii) above. 

2. Scattering of unbound stars 

2.1. Hyper-velocity stars 

Assume a MBHB of mass M = M\ + M 2 = 
M\{\ + q) (M 2 < Mi), and initial eccentricty 
e, embedded in a stellar cusp of mass Mb de- 
scribed by Maxwellian distribution function, 
with mean velocity cr. Quinlan (1996) showed 
that the hardening rate is relevant, and stays al- 
most constant, for separation smaller than 

GM 2 

known as the "hardening" radius of the 
MBHB. Here, it is important to notice that, 
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by measuring binary separations in units of a/,, 
scattering experiments results are independent 
on any pre-asigned value of <r. 

From the point of view of stars, repeated 
singshots results in a net heating of the distri- 
bution function, with a consequent erosion of 
the inner cusp. Though the overall energy bud- 
get is modest, as only a small fraction (< 1%) 
of the stars are in the loss cone at t = 0, nev- 
ertheless a substantial population of high ve- 
locity (often known as suprathermal) stars is 
produced. We define these stars by the condi- 
tion v > v esc , where the escape velocity v esc is 
taken as 

y esc = 2cr Vln(M B /M) * 5<r. (2) 

Equation |2 gives the escape velocity from the 
MBHB influence radius n = GM/lo 2 . The 
second equality comes from the adopted rela- 
tion M = 0.002 M B . 

An example of the effects of the slingshot 
mechanism on the stellar population is shown 
in Figure [2 where the initial (t = 0) velocity 
distribution (stellar mass per unit velocity) of 
interacting stars is compared to that after loss 
cone depletion, for an equal mass, circular bi- 
nary. It is clear that after the interaction with 
the binary, a large subset of kicked stars still 
lies in the (reduced) loss cone of the shrunk 
MBHB, i.e., they are potentially available for 
further interactions. We term stars that undergo 
secondary slingshots "returning". Note that a 
large fraction of returning stars have veloci- 
ties just below v esc . The fraction of secondary 
slingshots is relevant because the interaction 
with the MBHB, on one hand, increases the 
star velocity, while, on the other, moves the 
kicked stars on nearly radial orbits. 

Our calculations show that the high ve- 
locity tail of the distribution depends on the 
MBHB eccentricity, although the effect of 
changing e is small for small values of q. In 
this case, fewer stars are kicked out compared 
to the case q < 1, but, on average, at higher 
velocities. In general, both a small mass ratio 
and a high eccentricity increase the tail of high 
velocity stars. Integrating the curves over ve- 
locity gives the mass of interacting stars, which 
turns out to be 2M for equal mass binaries, 



1.2M for q = 1 /3 (e = 0.3), and =* 0.6M for 
q = 1/27 (e = 0.3). 

As already mentioned, not only the veloc- 
ity, but also the angular momentum of scattered 
stars is modified by the slingshot process. One 
of the more interesting effect is the flattening 
of the distribution of scattered stars into the or- 
bital plane of the binary. Moreover, scattered 
stars preferentially co-rotate with the binary. 
As a general trend, we can conclude that stars 
that acquire high velocities after the interac- 
tion, are preferentially moved into the binary 
orbital plane, on nearly radial, corotating or- 
bits. 

Recently, Levin (2005) pointed out that the 
anisotropy of the ejected star distribution is a 
decreasing function of time. We find a quali- 
tatively similar result in our scattering experi- 
ments, though the phyiscal context is quite dif- 
ferent. In figure|2 the time evolution of / 
for the ejected stars is shown, for different val- 
ues of q (assuming e = 0). As the binary or- 
bit decays, ejected stars are more and more 
isotropic. In fact, for small separations, V c » 
v esc , and even weak interactions can lead to fi- 
nal star velocities > v esc . The effect tends to 
be suppressed for small q, simply because for 
small mass ratios V c ~ v esc already at a — a/,. 

2.2. Mass ejection and coalescence 

Integrating the mass ejection rate J along the 
shrinking orbit, we can derive the ejected mass 
M e j. Note that, while J is indipendent on the 
time evolution of the orbit, M ej is an outcome 
of the hybrid model of the loss cone/orbit evo- 
lution. 

In Figure[3]results are given in terms of M e j 
normalized to M (left scale) and to Mo (right 
scale), as a function of q. Our results show 
that M ej ~ 0.7//, i.e., M ej /M ~ 0.7^/(1 + qf 
and M ej /M 2 ~ 0.7/(1 + q). Note that M ej /M 
(and M^jMi) is independent on the absolute 
value of M. We also checked that the value of 
the eccentricity is basically irrelevant for what 
concerns M e j. The bulk of the mass is ejected 
in a initial burst, though the effect is reduced 
for small q. The burst is due to the ejection of 
those stars present in the geometrical loss cone 
when the binary first becomes hard. For small 
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Fig. 1. Stellar velocity distribution (normalized to the total mass of interacting stars in units of the total 
binary mass M) for an equal mass, circular binary, at different stages of the hardening. The vertical line 
marks v esc a 5cr. Dashed lines: from top to bottom, distribution of stars in the shrinking loss cone before 
the 1st, 2nd, 3rd and 4th iteration. For clarity, the initial loss cone distribution is marked with a thicker line. 
Thin solid lines: from top to bottom, distribution of stars which have received 1, 2, 3 and 4 kicks. Thick 
solid line: final stellar velocity distribution after loss cone depletion is completed. 



q, mass ejection is already relevant at a a/,, 
as in this case V c is always > v esc . 

We may ask now if the amount of ejected 
mass is sufficient to shrink the MBHB orbit 
down to the GW-dominated regime. The time 
needed to reach coalescence emitting gravita- 
tional waves for a binary of mass M, initial ec- 



centricity <?, and separation a, is computed as 
in Peters (1964). 

Given that, it is straightforward to esti- 
mate the mass that needs to be ejected to have 
fcw(fl/) = 1 Gyr, i.e., the mass ejection needed 
to reach a final orbital separation at which GW 
emission leads to coalescence within 1 Gyr. 
Results are shown in Figure|5] The upper, dark 
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Fig. 2. Time evolution of angular momentum of ejected stars, in terms of L^/L^, for a circular binary, and, 
from bottom to top, for q =1/27, 1/9, 1/3 and 1 



shaded area defines such mass (in units of M) 
for a 1O 6 M MBHB. The upper area limit is 
in case of a circular binary, the lower limit for 
e = 0.9. The lower, light shaded area defines 
the same "mass that needs to be ejected" for a 
1O 9 M MBHB. Again, upper and lower lim- 
its are in the case e — and e = 0.9, re- 
spectively. Incidentally, the e = 0.9 limit of 
the light MBHB practically coincides with the 
e — limit of the heavy one. Comparing these 
estimates with our M e j/M curve (which, we 



recall, does not depend upon M nor e), we 
can conclude that, if stellar slingshot is the 
only mechanism driving MBHBs to aj, circu- 
lar binaries are not going to coalesce within 1 
Gyr after they first become hard. In the case 
of M — 10 6 M o , only extreme unequal mass, 
highly eccentric MBHB will be led to coa- 
lescence by GW emission (and hence, would 
be detectable by LISA). Higher masses are fa- 
vored, but, in any case, still relevant eccen- 
tricities and/or small mass ratios are needed. 



F. Haardt: Black Hole Binaries 



287 



1 - 



3 s 0.1 



0.01 



- 10 




Mg/M, 



Fig. 3. Ejected stellar mass M c j normalized to the total binary mass M (left scale, solid points), and to the 
mass of the lighter binary member M 2 (right scale, empty points), as a function of the binary mass ratio. 
The curves are polynomial interpolations. Note that M e j/M and M^/Mi does not depend on the absolute 
value of M, and are nearly independent on e. The upper, dark shaded area defines the mass (normalized to 
M) a M = 10 6 M binary needs to eject to reach a final separation 0/ such that few = 1 Gyr. Upper and 
lower limits to the area are for e = and e = 0.9, respectively. The lower, light shaded area is analogous, 
but for a M = 10 9 M Q binary. Incidentally, the e = 0.9 limit for the light binary practically coincides with 
the e = limit of the heavy one. 



Figure @] shows the binaries (at two different 
reference redshifts) in the M\ - q plane which 
would be resolved by LISA with SNR > 5, and 
the binaries (in the same plane) which are go- 
ing to coalesce within 5 Gyrs after loss cone 
depletion is completed. The overlap of the two 



areas gives the MBHBs which are potential tar- 
gets for LISA. As expected, the eccentricity is 
a crucial parameter, and it is clear how, even 
considering only slingshot of unbound stars as 
orbital decay driver, very unequal mass bina- 
ries are ideal LISA targets. 
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Fig. 4. In the plane Mi - q, the vertical shaded area limits LISA potential targets with SNR > 5. The 
diagonal shaded area on the lower right corner marks binaries that are going to coalesce within 5 Gyrs after 
loss cone depletion. In each panel, the reference redshift and eccentricity of the MBHBs are labelled. 



3. Conclusions 

We have made a first attempt of modeling the 
time dependent dynamical evolution of MBH 
pairs in stellar backgrounds using an hybrid ap- 
proach, i.e., using results of three-body scatter- 
ing experiment to follow the time evolution of 
the stellar distribution. Despite several approx- 
imations and limitations, we have obtained a 
number of results we briefly summarize in the 
following. 



A net outcome of the MBHB-stellar inter- 
action is a significative mass ejection, that we 
found scaling as M e j O.lfi, nearly indepen- 
dent on the binary eccentricity. This implies 
that, during every major merger, a stellar mass 
in the range 0.5 - 1M is ejected from the bulge, 
a possible clue to explain the observed mass 
deficit in ellipticals (Haenhelt & Kauffmann 
2002; Ravindranath, Ho & Filippenko 2002; 
Volonteri et al. 2003; Graham 2004). The prop- 
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erties of the ejected stars are well defined. They 
typically corotate with the MBHB, and show a 
large degree of flatness, as they are, preferen- 
tially, ejected in the binary plane. The ejected 
sub-population, if detected in nearby galaxies, 
would be an unambiguous sign of a relatively 
recent major merger. Though the total ejected 
stellar mass does not depend on the exact value 
of the binary eccentricity, the velocity function 
of the population does. In fact, an highly eccen- 
tric binary tends to produce a more pronounced 
tail of hypervelocity stars. 

The relevant quantity to asses coalescence 
is the ratio df/dow, where the final separation 
df is the orbital separation after loss cone de- 
pletion. In general, aj is reached in a short 
time, t ~ 10 7 yrs. Though rapid, slingshot is, 
in general, not sufficient to drive binaries to the 
gravitational wave emission regime, unless the 
binary is very massive (M > 1O 8 M ), and/or 
very eccentric (e > 0.7), and/or very unequal 
mass (q < 0.01, in fact the gap - aqw 
is lower for unequal mass binaries, though 
the shrinking factor is small). Our conlusion 
is supported by recent N-body simulations 
(Makino & Funato 2004; Berczik et al. 2005). 
In terms of LISA binaries (10 4 - 1O 7 M ), it is 
then important to study in details the possible 
role of other mechanisms, such as non spheric- 
ity of the stellar bulge and presence of circum- 
nuclear gaseous disks. 
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